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Abstract. In contrast to the canonical case, microcanonical thermodynamic func- 
tions can show nonanalyticities also for finite systems. In this paper we contribute 
to the understanding of these nonanalyticities by working out the relation between 
nonanalyticities of the microcanonical entropy and its configurational counterpart. If 
the configurational microcanonical entropy ui^(v) has a nonanalyticity at v = v c , then 
the microcanonical entropy lon{^) has a nonanalyticity at the same value e = v c of its 
argument for any finite value of the number of degrees of freedom N. The presence 
of the kinetic energy weakens the nonanalyticities such that, if the configurational 
entropy is p times differentiable, the entropy is p + \_N/2\ -times differentiable. In the 
thermodynamic limit, however, the behaviour is very different: The nonanalyticities 
do not longer occur at the same values of the arguments, but the nonanalyticity of 
the microcanonical entropy is shifted to a larger energy. These results give a general 
explanation of the peculiar behaviour previously observed for the mean-field spherical 
model. With the hypercubic model we provide a further example illustrating our 
results. 



PACS numbers: 05.20.Gg, 05.20.-y, 05.70.Fh 
1. Introduction 



In equilibrium statistical mechanics, nonanalyticities of thermodynamic functions are 
associated with thermodynamic phase transitions: more precisely, one commonly defines 
a phase transition point as the value of an external parameter (like temperature or 
pressure) where some thermodynamic function is nonanalyticjj] Such an identification 
is satisfactory in the canonical ensemble: as originally suggested by Kramers PQ, 

| Departing slightly from the standard definition, we use the notion of analyticity in the sense of a real 
function being infinitely-many times differentiable. This property typically, but not always coincides 
with the standard definition via the existence of a Taylor series. 
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nonanalyticities of thermodynamic functions calculated in the canonical ensemble may 
show up only in the thermodynamic limit N — > oo, where N is the number of degrees 
of freedom [2]. Moreover, such nonanalyticities separate different phases, i.e., regions 
of the parameters where the collective properties of the system are different. When 
studying physical models, one usually finds that thermodynamic functions have only a 
small number of nonanalyticities. 

In the microcanonical ensemble, however, the situation is different. First, 
nonanalyticities of the microcanonical entropy may even be present at finite N. In 
principle, this fact should have been known for a long time, because even one-degree- 
of-freedom systems like a simple pendulum or a particle in a double-well potential do 
show nonanalyticities of the entropy. Still, it came as a surprise to many researchers in 
the field to see how frequently such nonanalyticities are encountered in many-particle 
systems [3], IH El E] . The behaviour of such nonanalyticities as a function of the number N 
of degrees of freedom is remarkable: their number may grow with iV even exponentially, 
and their "strength" generically decreases linearly with N. The latter means that the 
first n derivatives of the entropy are continuous, where n is of order N Hj. Since 
the usual thermodynamic quantities of interest are given by low-order derivatives of the 
entropy, these nonanalyticities of order N can be observed only for very small N from 
noisy data. In the thermodynamic limit most of these nonanalyticities disappear. Those 
(if any) that survive are typically associated to thermodynamic phase transitions and 
coincide with the canonical nonanalyticities if equivalence of statistical ensembles holds. 
Therefore the finite-iV nonanalyticities are not easily associated to any phenomenon 
that one would call a phase transition in the usual sense, and it seems inappropriate to 
define phase transitions in the microcanonical ensemble only based on the presence of 
nonanalyticities of the microcanonical entropy. 

The general relation between nonanalyticities of the microcanonical entropy 
and thermodynamic phase transitions, i.e., the question which of the finite-system 
nonanalyticities survive in the thermodynamic limit, is a subject of active research (see 
e.g. [HI HQ])- It remains an open problem, despite some recent advances where — under 
suitable conditions — the "flatness" of stationary points was shown to be relevant to 
whether its thermodynamic limit contribution is vanishing or not [H]E]. m the present 
paper we will investigate the effect of a kinetic energy term on these nonanalyticities. 
The general properties of microcanonical nonanalyticities have been studied so far only 
when no kinetic term is present and the Hamiltonian is identified with the interaction 
potential energy. Standard kinetic energy terms in the Hamiltonian, i.e., quadratic forms 
in the momenta, are known to yield only trivial contributions to thermodynamics in the 
canonical ensemble. In the microcanonical ensemble a standard kinetic energy term may 
have a more noticeable effect: for instance, it may restore equivalence between canonical 
and microcanonical ensembles when only partial equivalence holds in the absence of a 
kinetic energy [12J. But adding a kinetic energy term has also a remarkable effect 
on nonanalyticities of the microcanonical entropy: For the exactly solvable mean-field 
spherical model it has been shown recently that, in the presence of a standard kinetic 
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energy term, a nonanalyticity of the entropy which for any finite N is located at a fixed 
value of the energy per degree of freedom, jumps discontinuously to a different value of 
the energy in the thermodynamic limit [6]. 

In this paper we argue that such a behaviour of nonanalyticities in the thermo- 
dynamic limit is not a peculiarity of that model, but rather a general property of all 
the nonanalyticities that survive as N — > oo. Moreover, we discuss the behaviour 
of nonanalyticities at finite N, showing that the kinetic energy weakens them: If a 
nonanalyticity of order m is present in the configurational entropy, adding a kinetic 
energy term to the Hamiltonian will increase the order roughly to m + N/2, i.e., the 
first m + N/2 derivatives of the entropy will be continuous. 

The paper is organized as follows. After giving some definitions and fixing some 
notation, in section [2] the behaviour of microcanonical nonanalyticities is discussed in 
general. More specifically, in section 12.11 we recall the results on nonanalyticities of the 
configurational entropy, while section 12.21 is devoted to the effect of the kinetic energy. 
In section |3] we present a simple example to illustrate this behaviour. We will finish 
with some concluding remarks in section HI 



2. Microcanonical entropy and its nonanalyticities 

We consider classical Hamiltonian systems with N degrees of freedom, with Hamiltonian 
function 7i : t— > R of the form 

1 - 

n(p,q) = -J2pt + V (Q) (1) 

i=i 

with some potential energy V : T N \— > R. We denote by Tjy C H N the iV-dimensional 
configuration space and by its cotangent bundle, i.e., the phase space. We shall 
usually denote configurations as q = (?i,---,9jv) and phase space points as 

(p,q) = (pi, . ..,p N ,qx, . . .,q N ) E A N . 

The fundamental quantity of the microcanonical ensemble is the microcanonical 
entropy as a function of the energy (per degree of freedom) e, 

s N (e) = —\nu N (e), (2) 

where 

u N (e)= [ dpdq5[n(p,q)-N6] (3) 
JA N 

is the density of states and 5 denotes the Dirac distribution^. A related quantity is the 
configurational microcanonical entropy as a function of the potential energy (per degree 
of freedom) v, 

s c N (v) = — lnu c N (v), (4) 

§ We define all thermodynamic functions per degree of freedom, which accounts for the factor l/N in 
the definitions. The Boltzmann constant /cb is set to unity. 
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where 

u%(v) = [ dqS[V(q)-Nv] (5) 

Jt n 

is the configurational density of states. The conflgur ational entropy equals the entropy 
when the Hamiltonian just consists of a configuration-dependent potential energy, 
Ji = V. This is often the case when studying discrete spin systems where a definition 
of conjugate momenta is difficult. 

Alternatively, one can define the integrated density of states 

Q N (e)= f dpdqe[Ne-H(p,q)}, (6) 
Ja n 

where is the Heaviside step function, and the corresponding entropy function is 

= 1 In n N (e). (7) 
Again, the configurational counterparts can be defined, where 

n c N (v)= [ dqe[Nv-V(q)] (8) 

and the entropy is given by 

a c N (v) = ^lnQ%(v). (9) 
The density of states (J3j) and the integrated density of states (E]) are related by 

UN{e) = —^p (10) 

and an analogous relation holds for the configurational quantities Q C N and uj c n . Under 
suitable conditions on the Hamiltonian the difference between the entropies vanishes in 
the thermodynamic limit, 

= lim s N = lim a N , s^ = lim s c N = lim a c N (11) 

N^oo N^oo N— >oo N— >oo 

(see section 3.3.14 of [13] for details). Unless explicitly noted, the general properties of 
entropy functions considered in the following hold true for both definitions (j2J) and (J7]). 

As mentioned before, the microcanonical entropy is not necessarily an analytic 
function, neither for finite iV nor in the thermodynamic limit. We will point out in 
the following that nonanalyticities of the entropy originate from stationary points of the 
Hamiltonian [8J. A stationary point of a function / : M C H n i— > R is a point x c e M 
such that df(x c ) = 0, and the value f(x c ) is called a stationary value of /. When 
stationary points are non-degenerate, i.e., the Hessian matrix Hess/ of / is nonsingular 
at all stationary points x c , the function / is called a Morse function. In this case all 
the stationary points are isolated. The index j of the stationary point x c is the number 
of negative eigenvalues of the Hessian at x c . Minima and maxima are stationary points 
corresponding to j = and j = n, respectively. 

In the following we will assume the potential energy V(q) to be a Morse function 
unless explicitly stated otherwise. Conceptually, this is an insignificant restriction, 
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since Morse functions on some manifold M form an open dense subset in the space 
of smooth functions on M [14J, and are therefore generic. Hence, if the potential is not 
a Morse function, we can deform it into a Morse function by adding an arbitrarily small 
perturbation!]]] 

If the Hamiltonian T~C(p, q) is of standard form (TjQ), its stationary points are of the 

form 

(Pc,g c ) = (0,g c ), (12) 

where q c is a stationary point of the potential energy V(q). Hence the stationary value 
of the Hamiltonian coincides with the stationary value of the potential energy, i.e., 
e c = v c . Let us first recall some results on the relation between the nonanalyticities of 
the configurational entropy s c N and the stationary points q c of V. 

2.1. Nonanalyticities of the configurational entropy 

The Morse property of the potential V ensures that its stationary points are isolated, 
so that we can safely restrict our attention to a single stationary point q c . We denote 
by v c = V(q c )/N the corresponding stationary value per degree of freedom. It has been 
shown in [8] that in this situation the configurational density of states ([5]) can be written, 
to the leading order, as 



vl , , _ P {v - Vc ) + hN ^ V V J = + [(v- v c ) N ' 2 - s ] (13) 
Vldet [Hess v (g c )]| 

with some 5 > 0, where P is a polynomial of degree smaller than N/2 in v — v c , and 

( (-l) j ^x (N - 2)/2 Q(x) for j even, 

h NJ (x) = I {-lf j+1)/2 x {N - 2)/2 TX- 1 In |x| for N even, j odd, (14) 

[ (_l)(A f -j)/2(_ ;E )(^-2)/2 e( _ ;E) for iv> ■ Qdd 

If there are further stationary points, the configurational density of states is given by the 
sum of the contributions of each stationary point. A proof (of an even stronger result 
including higher order terms) is given in [8]. This result can be rephrased as follows: 

(i) Every stationary point q c of V gives rise to a nonanalyticity of the configurational 
entropy s c N (v) at the corresponding stationary value v — v c — V(q c )/N. 

(ii) The order of this nonanalyticity is [(A^ — 3)/2j , i.e., s c N (v) is precisely [(A^ — 3)/2j- 
times different iable at v = v c ]^ 

Since the integrated density of states Q C N is obtained from uj c n by integration, the 
nonanalyticity of the entropy following from definition (J7|) is slightly weaker. In this 
case, statement (ii) has to be replaced by: 

(ii') The configurational entropy o~ c N {v) is precisely [(N — 1)/2J -times different iable at 

v = v c . 

|| In case V is not a Morse function due to a continuous symmetry, Morse-Bott theory [15] should allow 
to carry over essentially all the results from the theory of standard Morse functions. 
If With L^J we denote the largest integer smaller than x. 
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2.2. The role of kinetic energy 

We have pointed out in equation l fl2l) that, for standard Hamiltonians of the form (CQ), 
if q c is a stationary point of the potential energy V(q) then (0, q c ) is a stationary point 
of Ti.(p, q) and vice versa. Hence the kinetic energy is zero at stationary points and 
T~L{Vc, <?c) = V(q c ) for all stationary points (p c , q c ) of 7i. As a consequence, for all finite 
N the nonanalyticities of the configurational entropy — which we have traced back to 
stationary points in the previous section — show up at the very same stationary values 
as those of the entropy. As we shall see below, the presence of a kinetic energy term has 
a twofold effect on these nonanalyticities. 

1. At any finite N, the order of the nonanalyticity is increased by the presence of a 
kinetic energy term. 

2. More surprisingly, in the thermodynamic limit those of the nonanalyticities which 
survive jump to a different value of the energy. 

Both these results follow from the fact that for Hamiltonians of the class (CD) the density 
of states can be written as a convolution [13]. Defining a kinetic density of states as 



we can write the density of states as 

POO 

u N (e)= d 7 u4(7Kr(£-7) = / d 7 ^(e - 7)^(7), (16) 







where the configurational density of states u c N is given by equation (JSJ). 

Since equations ( fl5l) and ( fl6l) hold in the same form also for the integrated densities 
of states Q, all the following results will be valid also for integrated densities of states 
and entropies a defined as in equation (J7|). We will now discuss the finite- N case and 
the thermodynamic limit separately in the following two subsections. 

2.2.1. The finite-N case. At any finite N, the effect of the kinetic energy term on 
the order of the nonanalyticities of the entropy can be computed explicitly from the 
convolution integral (fTBl) . Such a calculation is reported in |Appendix A[ and the only 



additional input used is that — in accordance with the results from section 12. H — the 
nonanalyticities of the configurational density of states uj c n are of algebraic type. 

Alternatively, the result of that calculation can be anticipated via an intuitive 
argument: We know that u 1 ^ is smooth, and we assume for simplicity that the 
configurational density of states is nonanalytic only at v = v c and analytic elsewhere. 
As long as e < v c , it is evident from the right-hand side of (fl6l) that the nonanalyticity 
of u c N at v c is never reached in the integration. Hence the convolution integral ffl6l) is the 
integral over the product of two analytic functions and yields an analytic function. As 
soon as e > v c , the nonanalyticity of uj c n is inside the range of integration and induces a 
nonanalyticity in the convolution integral. As a simple example consider the convolution 
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integral (TTB|) with the choices u) N {x) = x and oj c n (x) = Q(x — a) with a > 0. Performing 
the integration yields 

u N {e)= l -{e-a) 2 Q{e-a). (17) 

As expected from the above reasoning, this function inherits the nonanalyticity at e = a 
from uo c N . This is in agreement with our previous observation that the stationary values 
of the Hamiltonian coincide with those of the potential energy. In this example, uj c n is 
a 0-times differentiate function, whereas is 1-times differentiable. The calculation 
in Appendix A shows that in general adding a kinetic energy increases the order of the 
nonanalyticity. 

The effect of a standard kinetic energy term on the nonanalyticities of the entropy 



at finite N, as computed in |Appendix A[ can be summarized as follows: 

(i) If the configurational density of states uj^(v) is nonanalytic at v = v c , then the 
density of states uj^i^) and the entropy SAr(e) are nonanalytic at e — v c ; 

(ii) The density of states un(e) and the entropy sn(e) at e = v c are differentiable 
|_iV/2j-times more often than the configurational density of states uj c n (v). 

As far as the entropy is concerned, statements (i) and (ii) above hold if v c is in the 
interior of the domain of the entropy. Both statements hold also for the integrated 
densities of states £In(e) and fl c N (e) and for the entropy cr N (e). 

2.2.2. Thermodynamic limit. We will assume in the following that the thermodynamic 
limit of the configurational entropy s c N exists, i.e., that u c N = exp(Ns c N ) increases 
exponentially with N asymptotically for large N. From definition (115j) it follows that 
the kinetic density of states ^(7) is related to the volume of an (N — l)-dimensional 
sphere with radius ^2N r y, 



N/2-1 



with 

•» = 2 fS)< 2iv)W/2 " (i9) 

As a consequence, this quantity likewise grows exponentially in iV and the thermody- 
namic limit of the kinetic entropy s% = \n(u^)/N exist. Under these conditions, the 
thermodynamic limit is known to exist also for the entropy s^ (see section 3.4.1 of [T3] ) 
and is given, apart from irrelevant constants, by 

1 X 

Soo (e) = lim — lnwjv(e) = lim — In max [0^(7) u c N {e - 7)] . (20) 

N— >oo iv N— >oo iv 7^0 

Denoting by 7(e) the value of 7 that realizes the extremum in the right-hand side of 
equation (1201 . this expression can be rewritten as 

Soo(e) = £[7(6)] + s^e - l(e)}- (21) 
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From the above equation it is apparent that 7(e) must be equal to the microcanonical 
average of the kinetic energy per degree of freedom, 

1 /I - \ 1 

^) = aK 2 = £ ~ N {V{q)) = £ ~ {V){£) > (22) 

which allows us to rewrite equation (1211) in the form 

*o«,(e) = a S } [e-(t;>( e )] + a So[(t;>(e)]. (23) 

If ^(u) is nonanalytic at v = t> c , then Soo(e) will be nonanalytic at e = £*, where e* is 
defined implicitly by 

(v)(e*)=v c . (24) 

Apparently e* differs from the value of v c , unless the average kinetic energy ( 1221) 
vanishes at e*. Hence, despite their common origin from the nonanalyticity of s c at 
v c , nonanalyticities of Soo(e) jump from their finite- N value of e to a different value in 
the thermodynamic limit. 

These properties of the nonanalyticities of the microcanonical entropy at finite 
and infinite N had previously been studied for a simple model system, the mean-field 
spherical model. The results reported in Refs. [6j [12] — surprising at the time — are all 
in agreement with and satisfactorily explained by the results reported in the present 
article. To further illustrate the predictions, we will discuss in the next section an even 
simpler, analytically solvable model which is of pedagogical value. 



3. A simple example: the hypercubic model 

The hypercubic model, introduced in [IB] , can be seen as an A-dimensional gener- 
alization of a particle in a one-dimensional potential made up of two square wells 
separated by a finite barrier. The Hamiltonian is of standard form (CQ), i.e., is given 
by a standard kinetic term plus a potential energy V(q). To define the potential, we 
consider a hypercube B of side length b, centered at the origin, as well as two disjoint 
hypercubes A + , A~ C B of side length a ^ 6/2, symmetrically arranged with respect 
to the hyperplane qi + qi + • • • + = (see sketch in Fig. [I]). The potential of the 
hypercubic model is then defined as 

{0 for qE {A + UA-}, 

Nv c for q E B\{A+ U A'}, (25) 
00 for q G H N \ B. 

The integrated density of states of this model can be computed analytically for arbitrary 
N, and the effect of a kinetic energy term will nicely illustrate the general results of the 
previous section. 
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Figure 1. Arrangement of hypercubes defining the potential energy (j^5T) of the 
hypercubic model for N = 2. 



3.1. Finite-N integrated microcanonical entropy 

The integrated density of states is calculated for the hypercubic model by inserting the 
potential (125'j) into definition 

N 



Jr" Ja+ua- V ~i 1 J 



+ 



dp 



dq<S) 



B\{A+UA~} 



N 



i=i 



r / N 2 \ r r N 

2a N / dpe[Ne- V ^- ) + (b N - 2a N ) / d P 6 N(s - v c ) - V 
•/ r-^ V ~^ 2 y y R iv [ ~^ 



AT r, 

Pi 



The remaining integrations of the form 



R^ 



2 N 2 



C7V7 



A' 



with C 



7T 



7V/2 



A? 



T(N/2 + 1) 



(26) 



(27) 



yield the volumes of iV-dimensional balls of radii 7 = \/2Ns and 7 = ^/2N(e — v c ), 
respectively. Inserting this formula into and making use of ([7]), the integrated 
microcanonical entropy 



a N [e 



N 



1 



1 



— + - ln2 + lna+ — \nC N + -(IniV + lne) 



N 



for < e < v r 



^ In C N + In [2a iV (2AA £ ) JV / 2 + (6* - 2a") [2N(e - v c )} N/2 ] for £ > v c , 
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can be computed. For the integrated configurational microcanonical entropy, a similar 
calculation yields 

{In a + — In 2 for < v < v c , 
N (29) 
In b for v ^ v c . 

From these equations we see that o~n and <j c N each have one nonanalyticity in the 
interior of their domains [0, oo), located at argument v c in both cases. The integrated 
configurational microcanonical entropy a c N is a piecewise constant function with a 
discontinuity^ at v — v c , whereas the functional form of the integrated microcanonical 
entropy o~n in the vicinity of v c is \s — f c | , in agreement with the general results of 
section 12.2.11 

3.2. Thermodynamic limit 

From the entropies (128]) and (1291) . the corresponding thermodynamic limit values 

0"oo(£) = hm <Jn{s) and ^(e) = lim a c N (e) (30) 
can be computed. Since 

lim — lnC N = - (1 + lnvr + ln2 - In iV) , (31) 

N-^oo N 2 

for < e < v c we find 

11 1 
<7oo(e) = lna + -lne + -lnvr + ln2 + -, (32) 

whereas for the case e ^ v c we can write 



t7 00 (e) = -(l + ln7r + ln2-lnJV) (33) 

+ 1 In {2a N (2Ne) N / 2 + {b N - 2a N ) [{2N) N / 2 {e - v c f/ 2 ] } . 

The argument of the logarithm in equation (1331) is the sum of two terms. Both terms 
are exponentially large in N, so that in the limit iV — ► oo only the larger one survives. 
There is a value 



<34) 

of e for which the two terms are equal, and apparently e* > v c . For e < e*, the first 
term in the argument of the logarithm survives, yielding the same functional form for 
o"oo as in the case < e < v c . For e ^ e*, however, the second term wins, and we obtain 
as a final result 

{- + \n2 + \na+ -Inn + -\iae for < e < £*, 

2 2 2 (35) 

- + In2 + ln6+-ln7r + - ln(e - v c ) for e ^ s*. 

+ Since the potential ((25]) is not a Morse function, it is no surprise that the nonanalyticity is not of 
the generic type given in equations (|13p and (TTJ 
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This expression shows that the entropy o' OQ (e) is analytic at e = v c but nonanalytic at 
e = e* 7^ v c . Therefore, although the entropy is nonanalytic at e — v c for all finite 
N, the nonanalyticity jumps to a different energy value e* in the thermodynamic limit. 
Furthermore, at this value the statistical average (v) of the potential energy per degree 
of freedom v equals the value of the finite- N nonanalyticity, 

(v)(e*)=v c . (36) 

A plot of the entropy for various finite values of N as well as in the thermodynamic 
limit is shown in Fig. [2j Again, the results for the hypercubic model confirm the general 
reasoning of section 12.2.21 



cr N 
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5 



Figure 2. Entropy ctjv of the hypercubic model, plotted as a function of the energy e 
for N = 2 (red), N = 8 (green), N = 32 (blue) and N = oo (black). The numerical 
values of the parameters are v c = 1, a = 1, b = 2, so that e* = 4/3. 



From equation (I35J) the microcanonical temperature 

/ds\r\ - l { 2£ for e < £*, 

Tie) = (—Z) = \ ' (37) 

V ; \ de J \ 2(e - v c ) for e > e*. V ; 

is easily computed. The jump in temperature at e* is a signal of a microcanonical 
first-order phase transition. 

4. Concluding remarks 

Previously, nonanalyticities of the microcanonical density of states and of the microca- 
nonical entropy usually had been investigated for configurational quantities. This choice 
was mostly motivated by the belief that the effect of a standard kinetic energy term, 
i.e., a quadratic form in the momenta, is trivial anyway. However, we have shown in the 
present article that the effect of a standard kinetic energy term on the nonanalyticities 
of the microcanonical entropy is quite remarkable, both for finite and for infinite N. 
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If the configurational microcanonical entropy oj c n (v) has a nonanalyticity at v — v c , 
then the microcanonical entropy ujn(e) has a nonanalyticity at the same value e — v c of 
its argument for any finite value of the number of degrees of freedom N. The presence 
of the kinetic energy weakens the nonanalyticities. More precisely, if the configurational 
entropy is p times differentiable, then the entropy is p + L/V/2j-times differentiable. In 
the thermodynamic limit, however, the behaviour is very different: The nonanalyticities 
do not longer occur at the same values of the arguments, but the nonanalyticity of the 
microcanonical entropy is shifted to a larger energy as given by equation ( pMl) . 

These results give a general explanation of what had previously been observed for 
the mean- field spherical model [6]. With the hypercubic model we have provided a 
further example for which both, the configurational microcanonical entropy and the 
microcanonical entropy can be computed analytically. Due to its simplicity, this model 
is of pedagogical value and nicely illustrates the general behaviour of nonanalyticities 
of the entropy. 

The study of nonanalyticities of the configurational microcanonical entropy, their 
relation to stationary points of the potential energy, and their connection with topology 
changes of the constant-potential surfaces in configuration space, has attracted some 
interest recently and has proved useful to obtain a deeper understanding of phase 
transitions (see [Him [10] f° r reviews). With the present work, we add another piece to 
this understanding by providing the relevant "translation rules" between configurational 
quantities and their total-energy counterparts. 
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Appendix A. Finite-iV nonanalyticities of the density of states: explicit 
calculation 

We now explicitly calculate the convolution integral ( 1161) at finite N. For the sake of 
simplicity, let the configurational density of states oj%{v) have just one nonanalyticity 
at v — (this choice of the stationary value will ease the notation). According to the 
general results reported in [8] and recalled in section [21 we assume, apart from irrelevant 
multiplicative constants, 



where the distribution forces u c N to vanish when the potential energy per degree of 
freedom v is smaller than its minimumE] v min , is smooth and u±(v) = O (v p ) is 

* Since we assume v c = 0, we have v m i n < 0. 




(A.l) 
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the nonanalytic part; more precisely, to the leading order 

for v < 
for v > 





\v\- 




'\v\- 



_ , p for v < 0, 

-'±<''> : { -l L ,. (A.2) 

where c~,c + G R and 2p G N . With the the explicit expression ( fl8|) for the kinetic 
density of states the convolution integral tTTHj) becomes 

/>oo 

u N (e) = u N (e) + a N d-y <~i N,2 ~ l u)±(e - 7) (e - 7 - v min ) (A. 3) 

./o 

with some smooth function u)n(e). We are interested in the integral containing the 
nonanalyticity, 



d 7 7 AV2 -W(£ - 7) © (e - 7 - ^min) , (A.4) 
and we can write I = I\ + I2 with 

/■oo 

h= d7 7 Ar/2 - 1 c + |£-7| p 6(£-7), 



(A.5) 



00 



d7 7 AV2 ~ 1 c~k - 7l P e( 7 - e) 6 (e - 7 - iw) . (A.6) 
With the substitution x = j/e we obtain 

jj = c + £ P +Ar / 2 0(e) / dxx^-^l - x) p = c + £ p+JV/2 fi(Ar/2, 1 + p) 6(e), (A.7) 

where B(u,w) = T(u)T(w)/T(u + w) is the beta function. 

For the calculation of I2 it is convenient to treat the cases e < and e > separately. 
For e > 0, 

d7 yv/2-i (T _ £)P) (A8) 

and with the change of variables x = e — 7 one gets 

lin 

dx {-xf (e - x)^ 1 . (A.9) 
The substitution y = x/e yields 



/ 2 = _ c -(_iyV>+^/2 



dy^l-y)^ 8 - 1 

= -c-(-l)%^/ 2 S We (l + p,iV/2), (A.10) 

where S z denotes the incomplete beta function Jj] For e < 0, a similar calculation gives 

h = c^-ir^i-eY^B^j^N^, 1 + p)6(e - tyj. (A.ll) 

Assembling the pieces together we finally obtain for the nonanalytic part of the density 
of states 

u N (e) - u> N {e) oc \e\ p+N/2 [b~{e) 0(-e) + 6+(e) 6(e)] , (A.12) 

jj The incomplete beta function B z has a branch cut discontinuity in the complex plane running along 
the negative real axis. This implies that, when expressing I2 in terms of B z , the correct branch has to 
be used. The given integral representations, however, are unambiguous. 
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where 

6"(e) = c-(-l)- Ar / 2 5 1 _, ;min/£ (iV/2, 1 + p)8(e - *w), (A.13) 
6+(e) = c + B (N/2, 1 +p) - C -(-l)^ mm/e (l + p, N/2). (A.14) 

The function given by equation (1A.12I) . and thus the entropy Sn(e), is nonanalytic at 
e = = v c . At this point, sn is [p + (N — 1)/2J -times differentiable, whereas the 
configurational entropy s c N we started out with is only \_p — 1/2 J -times differentiable. 
Without any modifications the calculation can be carried over to the integrated density 
of states Qjy. 
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